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We calculate the cross section for the exclusive production of J PC — ++ glueballs Go in asso- 
ciation with the J/tp in e + e~ annihilation using the pQCD factorization formalism. The required 
long-distance matrix element for the glueball is bounded by CUSB data from a search for resonances 
in radiative T decay. The cross section for e + e~ — > J/ip + Go at y/s = 10.6 GeV is similar to ex- 
clusive charmonium-pair production e + e~ — > J/ip + h for h = n c and XcO, and is larger by a factor 
2 than that for h — rjc{2S). As the subprocesses 7* — » (cc)(cc) and 7* — > (ccj(gg) are of the same 
nominal order in perturbative QCD, it is possible that some portion of the anomalously large signal 
observed by Belle in e + e~ — > J/tpX may actually be due to the production of charmonium-glueball 
J/4<Gj pairs. 



Bound states of gluons provide an explicit signature 
of the non-Abelian interactions of quantum chromody- 
namics. In fact, in a model universe without quarks, 
the hadronic spectrum of QCD would consist solely of 
color-singlet glueball states. In the physical world, the 
purely gluonic components mix with qq pairs, leading to 
an enriched spectrum of isospin-zero states as well as qqg 
hybrids. The existence of this exotic spectrum is as es- 
sential a prediction of QCD as the Higgs particle is for 
the electroweak theory. 

Lattice gauge theory predicts the spectrum and quan- 
tum numbers of gluonic states. According to a recent cal- 
culation [l|, the ground-state masses for the J PC = ++ 
and 2 ++ glueballs Qj are 1.73 and 2.40 GeV, respec- 
tively.Thus far, the empirical evidence for glueballs is not 
decisive, probably because of complications from mixing 
with the quark degrees of freedom, but there are indica- 
tions of an extra neutral scalar state perhaps due to a 
glueball of mass (before mixing) near 1.7 GeV 0. 

An important mechanism for producing glueballs is 
the radiative decay of heavy quarkonium, particularly 
J/ip — > jQj and T — > jGj In these reactions, the 
quarkonium decays to an intermediate jgg state which 
then can couple to any charge conjugation parity C = + 
isospin 7 = gluonic or hybrid state. For example, 
the BES Collaboration |4J has observed the radiative de- 
cays of the J/ip and the ip(2S) to 7/o(1710), a glueball 
candidate. In this Letter we focus on another optimal 
mechanism for the production of Qq and Q2 at e + e~ col- 
liders, the reaction e + e~ — > 7* — * HQj, H = J/ip or 
T 0, in which a C = + glueball can be produced in 
association with a quarkonium state from the subpro- 
cess 7* — > (QQ)(gg)- Two-gluon components in 77 par- 
ticles have been estimated recently [6]. One of the six 
Feynman diagrams for the subprocess is shown in Fig.Q 
the remaining diagrams are permutations of the photon 
and the two gluons. A related reaction 7* — > ir°Qj has 
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FIG. 1: Feynman diagram for 7* — > H + Gj- 

been considered Q as a source of pseudoscalar glueballs. 
We shall show that these reactions satisfy perturbative 
QCD (pQCD) factorization. Unlike radiative quarko- 
nium decay, this channel imposes no a priori limit on 
the mass of the glueball. 

The main background to charmonium-glueball pro- 
duction e + e~ — > J/tpQj is exclusive quarkonium pairs 
such as 7* — * J/iprjc, arising from the subprocess 7* — > 
(cc)(cc). The exclusive production of charmonium pairs 
has in fact been observed recently with a substantial 
rate at Belle 0. The rates for exclusive charmonium- 
pair production reported by Belle are sig nificantly larger 
than predictions based on pQCD [SUlfJ. The Belle ex- 
periment identifies one member of the pair, the J/ip, 
via its leptonic decay; the other quarkonium state is 
inferred by identifying the missing mass of the spec- 
tator system with the charmonium states rj c , Xco, an d 
?7c(25) which occur within the detector mass resolution. 
As noted in Refs. [ill ll2T ]. some of the Belle signal for 
quarkonium pairs may be due to two-photon annihilation 
e + e~ — > 7*7* — > J/tpJ/tp. Here we note that because 
the subprocesses 7* — > (cc)(cc) and 7* — > (ccj(gg) are of 
the same nominal order in pQCD,it is possible that some 
portion of the signal observed by Belle in e + e~ — > J/ipX 
may actually be due to the production of J/ipQj pairs. 
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In general, exclusive amplitudes can be computed 
in QCD by convoluting the light-front wavefunctions 
ipn/H { x ii kj_i) of each hadron with the correspond- 
ing n-particle irreducible quark-gluon matrix elements, 
summed over n |l3j . For hadronic amplitudes involv- 
ing a hard momentum transfer Q, it is usually possi- 
ble to expand the quark-gluon scattering amplitude as a 
function of h^/Q 2 . The leading- twist contribution can 
then be computed from a hard-scattering amplitude Tu 
where the external quarks and gluons associated with 
each hadron are collinear. Furthermore, only the mini- 
mum number of quark and gluon quanta contribute at 
leading order in 1/Q 2 . In our case, the relevant hard- 
scattering amplitude is Th("i* — > ccgg) computed with 
collinear c and c and collinear gg. As Th at leading twist 
is independent of the constituent's relative transverse mo- 
mentum kj^i, the convolution with the light-front wave- 
functions and the integration over the relative transverse 
momentum then projects out the L z = component 
of the light-front wavefunctions with minimal n - the 
hadron distribution amplitudes 4>h(x,Q). 

In this Letter we shall calculate the cross section for 
e + e~ —* HQj = o.2 using pQCD factorization. The ampli- 
tude at leading twist can be expressed as a factorized 
product of the perturbative hard-scattering amplitude 
Th{i* — ► QQgg) convoluted with the nonperturbative 
distribution amplitudes for the heavy quarkonium and 
glueball states. We shall find that 7* — * J/tpGo produc- 
tion dominates over that of J/1PQ2, and show how the 
angular distribution of the final state can be used to de- 
termine the angular momentum J and projection J z of 
the glueball. We shall show that only J z = ±2 tensor 
states are produced by the pQCD mechanism at leading 
twist. A bound on the normalization of the distribution 
amplitude for the glueball state can be extracted from a 
resonance search by CUSB in T — * 14]. We shall 
show that the rate for e + e~ — > J/tpGo production could 
be comparable to the corresponding NRQCD prediction 
for e + e~ — * J/tprjc without exceeding the CUSB bound 
from radiative T decay. 

The distribution amplitude 4>h(x,Q) required for the 
formation of the H in a hard process is directly related 
to the NRQCD matrix element for the leptonic decay 
rate of H. Its x dependence is peaked at x ~ 1/2. The 
key quantity which determines the normalization of the 
7* — > HQj processes is then the distribution amplitude 
(f>.j(x, Q) of the Qj. The pQCD factorization picture pro- 
vides a direct relation among the various glueball pro- 
duction processes, as they all involve the same process- 
independent (f>j(x,Q). The (f>j(x,Q) can be determined 
phenomenologically by fitting to the measured produc- 
tion rate of a glueball candidate. In leading- twist ap- 
proximation the spin structure of the two-gluon system 
in hard-scattering amplitude becomes that of a massless 
spin- J = 0, 2 state. Therefore the field-theoretic defini- 



tion of the 4>j(x, Q) in light-cone gauge reduces to [If 



4>j(x,Q) = 



F J 
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^2{N 2 -I) J {2nfk+x{l-x) 

x e ~i(xk + z~ -kj_-zj_) 

x(gj\T G+°*(0+,z-,z x )G+P(0)\0), (1) 



where x and are the light-cone momentum fraction 
and transverse momentum of a gluon inside the Qj with 
momentum k = (k + = n ■ k,k~ = h ■ k,0±) and mass 
Mgj. The S— wave component is projected out by in- 
tegrating over kj_. The light-like vectors n and n sat- 
isfy n 2 — n 2 = and n ■ n = 2. The tensor 
projects the massless spin-J components; they are de- 
fined by F°p = [-g a p + \ {n a fif3 + npn a )]/^/2 and F 2 is 
the massless spin-2 polarization tensor e a p . The glueball 
distribution amplitude can also be defined from the two- 
gluon light-front wavefunctions ipgj (x, k^, Aj) with gluon 
spin projection Ai = Sf = ±1, integrated over transverse 
momentum in light-cone gauge A + = 0. 



The relative rates for the production of heavy scalar 
glueballs with higher radial number N are determined by 
the normalization of the corresponding glueball distribu- 
tion amplitudes. In effect, the integral of the distribution 
amplitude over x is the relativistic generalization of the 
Schrodinger wavefunction at the origin. Thus the dis- 
tribution amplitudes for the ++ glueballs tend to scale 
inversely with their mean radius (rjv)- According to bag 
models (tm) ~ 0.6 fm, independent of N, suggesting 
equal rates for the heavier glueballs. On the other hand, 
the virial theorem extended to the light-front formalism 
suggests that mean transverse momentum and l/(rjy) 
increase monotonically with glueball mass. If this is the 
case, then the production rate in the 7* — > HQj will tend 
to increase for heavier glueball states, assuming that the 
annihilation energy ^/s poses no phase-space restriction. 
Lattice gauge theory and light-front Hamiltonian meth- 
ods should eventually determine the glueball distribution 
amplitudes, thus providing consistency checks on the pro- 
duction mechanisms considered here. 

As noted above, the amplitude for 7* — * H(p)Qj(k) 
can be computed as the convolution of Th{"i* — > QQgg) 
with 4>j(x,Q) weighted by the NRQCD matrix element. 
In leading twist fc~ is neglected and, thus the glueball 
momentum is approximated by k = k + n/2 in Th(j* — > 
QQgg). The resulting effective vertex Aj is 0] 
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where [i and e# are the vector indices for the 7* and 
polarization vector for the H, respectively. The mass 
and fractional charge of the heavy quark Q are ex- 
pressed as m,Q and eQ. Here {Oi)h and (Q\)h are 
the vacuum-saturated analogs of NRQCD matrix ele- 
ments (Oi( 3 Si)}h and (Q]( 3 S-t))rr for annihilation de- 
cays defined in Refs. and 0], respectively. To lead- 
ing order in the heavy-quark velocity vq in the quarko- 
nium rest frame, the (Oi( 3 Si))h is related to the ra- 
dial wavefunction at the origin i?(0) in the color-singlet 
model [lflj l and the decay constant which is defined 
by (0\J£ m .\H) - 2M H e Q f H e^, as {O x ) H = f#|i?(0)| 2 = 
2MjjJh- The nonperturbative factors for Gj are written 
as I = f Q dx<j) (x, Q) and J 2 = J* Q dx4> 2 (x, Q)/[x(l — x)]. 
In leading twist the valence gluons are collinear and 
therefore the only allowed polarization states for Gi are 
A2 = ±2. For Gi production the longitudinal polariza- 
tion is prohibited by Bose symmetry. This is true for any 
production process for G2, for which pQCD factorization 
is valid. The amplitude © for Q2 is proportional to the 
factor (<2i)fj/mQ which is suppressed to (Oi)#/mg by 
Vq. Therefore, in the remainder of this Letter we only 
consider Go', the analysis for Q2 can be found in our forth- 
coming publication 15]. Using the vertex @, we obtain 
the width Tq for radiative T decay into Go as 



r„ 



where $q 
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In Ref. |20j, the decay rate 



for the process T — > 7/0 has been calculated treating 
fo as a glueball candidate. The r agrees with Eq. (5) 
of Ref. |20|, after including a missing factor 2/3 and ne- 
glecting m Go mi. 



Our result for the differential cross section for 



J/ipGo normalized to 



47ra 2 /(3s) 



i.e., aligned with the lepton beams, then only final states 
with J z = ±1 can contribute, because the e + and e~ 
annihilate with opposite chirality. Thus in the case of 
scalar glueballs, the J/ip with helicity ±1 is produced 
with a 1 + cos 2 6* distribution. If the J/ip is longitu- 
dinally polarized, the cross section must vanish in the 
forward direction, and thus it has a sin 2 8* distribution. 
The rate integrated over angle is 



R 



J/ipGo 
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The size of the cross section can be estimated using the 
asymptotic form of the ratio 1Z = Rj/jpg /Rj/ t f, nc 



n 



m c |Jo| 



r 2 (l-r 2 )(l-^) 2 <Oi), c 



(8) 



is 



where we neglected QED contributions to J?j/^ r)c given 
in Ref. In the ratio 1Z the phase-space factor $q 6 
cancels the yl — r 2 for e + e~ — * J/iprj c . The a s 's for the 
two processes are written distinctively because they have 
different effective scales. However, the main uncertainties 
from the choice of running coupling scale and scheme 
largely cancel in the ratio TZ. Here ™q[° represents the 
ratio of the square of the wavefunction at the origin of 
the glueball compared to that of the rj c . 

We next investigate whether some portion of the 
anomalously large signal for J/ip + rj c , XcO, and r/ c (2S) 
observed by the Belle Collaboration could actually be 
coming from the process e + e~ — * J/ipGo- We calculate 
the cross section assuming glueball mass Mg the same 
as those for rj Cl Xco, and r] c (2S). In order to predict the 
production cross section <rj/ipg , we need to know the 
nonperturbative factors (Oi)jm and Iq. The {0\)j/^, 
is determined through the leptonic decay rate of J/ip- 
As the glueball distribution amplitude is process inde- 
pendent, we can extract an upper bound to Iq from the 
CUSB data for the resonance search from T — > jX. We 
follow the method used in Ref. [13] • The branching frac- 
tion Br^CJo] for the process T — ► jGo is obtained by 
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where 9* is the scattering angle in the center-of-mass 
frame, r = 4m c /y/s, and the phase-space factor <&3 e is 
defined by 



- {Mj N , + Mg ) 2 ] [s - (M j/ip - Mg ) 2 ] .(6) 



The angular factors in the expression JSJ can be under- 
stood physically. If the hadron pair is produced at 9* = 0; 



Br[ 7 £ ] = 



r[e+e ] NRQCD 



x Br[e + e ] c 



(9) 



where Br[e + e~] oxp . = 2.38% and r[e^e — ]nr,qcd = 
27re 2 a 2 (Oi)x/(3m 2 ). In the ratio (Oi)r dependence 
cancels. The branching fraction must be less than al- 
lowed by the CUSB excluded region. In order to extract 
the bound, we note that the mass resolution of the CUSB 
data is 20 MeV. If the decay width T[Go] of the Go is larger 
than the resolution, one must rescale the boundary of 
the excluded region by the factor r[£ ]/20 MeV. The de- 
cay width T[Go] cannot be computed using perturbation 
theory because factorization is not valid for this non- 
perturbative quantity. However, if Belle's J/ipr/ c signal 
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TABLE I: Upper limits to the nonperturbative constant [io| , 
cross section aj/^g , and the ratio u.j/^g /uj/^, h at y/s — 
10.6 GeV, assuming Mg = Mh, where h = r) c , XcO, and 
r) c (2S). The limits are determined by the T — > 7X search of 
the CUSB Collaboration Q. 



Mg Q = M h 


h = r]c 


XcO 


Vc(2S) 


|/o|Lk (10- 3 GeV 2 ) 


5.2 


5.8 


6.2 


max 


1.4 fb 


1.5 fb 


1.6 fb 


„max /_ 

cr j/^g / a J/i'h 


0.63 


0.72 


1.9 



also contains J/i/jQo, T[Qo] must be less than 110 MeV, 
which is the full width at half maximum of the r\ c peak 
in the Belle fit to the J momentum distribution. The 
first row in Table d gives the upper limits to |io| 2 for 
m b = 4.73 GeV and Mg = 2.98, 3.42, and 3.65 GeV 
corresponding to Mg = M Va , M Xc0 , and M v ^ 2 s), re- 
spectively. Values for |Jo| 2 above the bound are excluded 
by 90% confidence level. We choose a s (fi 2 ) = 0.26 us- 
ing the modified minimal subtraction (MS) scheme and 



the scale fi 



o-5/3 



1231 where ( k ) is the mean 



3-momentum squared for a single gluon. 

Now we are ready to find upper limits to <J,j/^,g at 
the B-factories. Substituting |/olmax to Eq. (JZJ), we get 
the cross sections <Jj/^g — R.j/ipg a u+u- in the second 
row in Tabled In order to make our prediction consistent 
with the previous analyses on exclusive charmonium-pair 
production we use the same input parameters given in 



Refs. 9, 



0103: (Oil 



J/i/> 



0.335 GeV J , m c = 1.40 GeV, 



and Mj/^ = 3.10 GeV. The strong coupling constant is 
chosen to be a s = 0.260, 0.264, and 0.265 for Mg = 
M Vc , M XoQ , and M^^S) ) respectively, applying the same 
method used for the radiative T decay. The ratio to the 
cross sections for exclusive charmonium-pair productions 
are given in the third row in Tabled 

The cross sections for J/ip + rj c , XcO, and r) c (2S) re- 
cently measured by the Belle Collaboration are not well 
understood within NRQCD. Based on the assumption 
that the measured signals at Belle include the J/ip + Go 
signal within the mass region corresponding to rj Cl XcO, 
and T] C (2S) we get the cross section for J/i/jGo- We thus 
find that the upper limit to the cross section <Jj/.^g is 
comparable to the NRQCD prediction of the cross sec- 
tions for e + e~ — > J/ip + h for h = r/ c and XcO, and larger 
by factor 2 to that for h — r] c (2S), suggesting the possi- 
bility that a significant fraction of the anomalously large 
cross section measured by Belle may be due to glueballs 
in association with J/ip production. In fact, there is a 
possibility of a resonance signal in the Belle data for 
e + e _ — > J/4>X at the missing mass M.x ~ 1.7 GeV. 
A resonance search in the radiative T(nS) decay by the 
CLEO Collaboration and an independent study by the 
BaBar Collaboration on charmonium-pair production in 



e + e annihilation will provide stringent tests of this sce- 
nario. 
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